Quantum discord and quantum phase transition in spin chains 
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Quantum phase transitions of the transverse Ising and antiferromagnetic XXZ spin S = 1/2 chains 
are studied using quantum discord. Quantum discord allows the measure of quantum correlations 
present in many-body quantum systems. It is shown that the amount of quantum correlations 
increases close to the critical points. The observations are in agreement with the information 
provided by the concurrence which measures the entanglement of the many-body system. 
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I. INTRODUCTION 

Quantum phase transitions (QPT) occur when the 
ground state of a many-body system at absolute zero 
temperature undergoes a qualitative change by varia- 
tion of a coupling and/or an external parameter 1 -. QPT 
are present in spin systems as for example in the one- 
dimensional XY model 2 -^ and in the antiferromagnetic 
XXZ model^^. Quantum phase transition in spin sys- 
tems have been widely characterized considering the pair- 
wise entanglement between two spin sitesi. The concur- 
rence which is a good measure of entanglement^ is max- 
imal close to the critical points and its derivatives can 
signal more precisely the presence of a quantum phase 
transition at the critical points^. 

In another register Ollivier and Zurek introduced the 
quantum discord as a measure of quantum correlations 
between two separated sub-systems of a many-body 
quantum system- -. In information theory the mutual 
information is defined as the measure of correlation be- 
tween two random variables A and B. For classical sys- 
tem the mutual information can be expressed in terms 
of two equivalent expressions T{A; B) and J{A\ B) (ex- 
plicit definitions of each expressions of the mutual in- 
formation are provided later in the text). However for 
quantum systems the two expressions are not equal and 
their difference gives the quantum discord. In parallel 
to the concurrence, quantum discord measures the quan- 
tum correlations of a quantum system. In contrary to 
the concurrence, the quantum discord can be different 
to zero even if a quantum systems is separable. The 
Werner state, pv/emer = ^-j^l + Al^)^! where A is a 
parameter and \^f) = (|00) + |11)) /v2) provide an exam- 
ple for which the concurrence is equal to zero (separable 
state) for A < 1/3 while the quantum discord is different 
to zeroi^. Hence quantum discord shows the existence 
of quantum correlations where concurrence shows that 
the system is separable. We point that the concurrence 
measures the nonlocal quantumness of correlations while 
quantum discord measures the total amount of quantum 
correlations of a state p. 

Quantum discord is used in several contex t 11 ' 12 for its 
property to discernate quantum and classical correlations 
of a quantum systems. The behaviour of the quantum as 



well as the classical correlations close to the critical points 
of QPTs can be studied. We focus our attention here 
on the quantum phase transitions of the one-dimensional 
Ising model which derives from the XY model and on 
the antiferromagnetic XXZ chain. Using the quantum 
discord we show that the amount of quantum correlations 
increases in the region close to the critical points. In 
opposite the amount of classical correlations close to the 
critical points decreases in the XXZ model while it is 
simply monotonous in the Ising model. 

The outline of the paper is the following. In section |TT] 
we recall the elements of information theory which leads 
to the construction of quantum discord as demonstrated 
by Ollivier and Zureki^. In section IIIII and IIVI the quan- 
tum correlations close to the quantum phase transition 
of the Ising and XXZ models are studied. Appendix [3] 
provides the derivation of the quantum discord for spin 
5=1/2 chains. 



II. QUANTUM DISCORD 

Information on the correlations between two random 
variables A and B can be obtained by measure of the 
mutual information. In classical information theory the 
mutual information reads 



I (A; B) = H(A) + H(B) - H(A, B), 



(1) 



where the information entropy H(A) = 
PA=a and PA=a is the probability 
that a is the realisation of the random variable 
A. The joint entropy of A and B is defined by 
H(A,B) = -J2a,bPA=a,B=b\ogp A =a,B=b where 
PA=a,B=b is the joint probability of the variables A 
and B to take respectively the values a and b. Using 
the Bayes rule PA,B=b = PB=b x PA\B=b the mutual 
information can be rewritten into the equivalent 
expression 



J(A;B) = H(A) - H(A\B), 



(2) 



where H(A\B) = -J2a,bPA=a\B=b^ogp A =a\B=b is the 
conditional entropy of the random variables A and B 
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and PA= a \B=b is the conditional probability of a being 
the realisation of the random variable A knowing that b 
is the realisation of the random variable B. For classi- 
cal random variables the two expressions |T]) and Q are 
strictly equivalent. However for quantum systems the 
two expressions of the mutual information are in gen- 
eral not equivalent and leads to define the quantum dis- 
cord 6 (A ■ B) = 1(A;B) — J{A;B) as demonstrated in 
refs i 10 ^ 1 . Quantum discord is the difference between the 
two classicaly equivalent expressions of the mutual infor- 
mation. 

The mutual information expressions (JTJ) and ^ for 
quantum systems are obtained by replacing the Shan- 
non entropy by the von Neumann entropy. The von 
Neumann entropy for the joint state reads H(A,B) = 
~^ T A,BPA,B m /°^,e where the trace runs over the Hilbert 
spaces of the sub-systems A and B, pa,B is the density 
matrix of the quantum system. In a similar way the 
entropy of the system A reads H(A) = — T^aPa^°EPA 
where pa is the reduced density matrix obtained by tak- 
ing the trace over all states of the system B, pa = 
TtbPA,B- The quantum conditional entropy H{A\B) in 
equation @ quantifies the ignorance of the state of A 
knowing the state of B. It is worked out from the density 
matrix Pa\ii b which is related to the density matrix of 
the system pa.b through 10 

PA\uf = KjPA^Rj/Pj, 

where U B is the projector onto the state j of the 

sub-system B and the probability pj = ^a,B^PA,B- 
The quantum conditional entropy then reads H (A\B) = 

tf(^|{nf}) = £ ift #(^ |nf ). 

Finaly the quantum discord measuring the quantum 
correlations between the two quantum sub-systems A 
and B reads 

8{A:B)= mm {h(A) - H(A,B) + H(A\ {llf })}. 

It must be noticed that the quantum discord is minimized 
over the set of state {rrj}. The information obtained 
on the sub-system A is affected by the measure on the 
state of the sub-system B. In order to get the maximum 
information on the sub-system A we need to consider 
the projection onto the states of B that disturb least the 
overall quantum system. Maximize the information on 
the sub-system A correspond to minimize the quantum 
discord with respects to the set of projector {nj}. 

Quantum discord provide information on the quantum 
nature of the correlations between two systems. If two 
sub-systems A and B are correlated classicaly the quan- 
tum discord is equal to zero. Moreover quantum discord 
shows that quantum correlations can be present for states 
that are not entangled. For example in Werner states the 
quantum correlations are still present while the system is 
separable^. To be more precise we can consider a quan- 
tum mixed state p = \ip)(tp\ acting in an Hilbert space 



Ha <8> Hb ■ The quantum mixed state is separable if its 
entanglement is equal to zero. In this case the quan- 
tum separable state can be written as a product of quan- 
tum states and reads p — Y^iPil^hAti^hAl ® \^i,b){^i,b\ 
where | V*,-4.) e %Ai H>i,B) € Hs and Y^iPi = 1 - Tne two 
systems A and B can be correlated through the separa- 
ble state p and the quantum states \ipi,A) an d IV^b) do 
not have in general any classical counterpart. A conse- 
quence of the quantumness of the states \iPa) and \4>b) is 
that the correlations between the two systems A and B 
have quantum as well as classical nature. The quantum 
discord is non-zero while the entanglement of the quan- 
tum separable state is equal to zero. A good example is 
furnished by the two-qubit separable mixed stated 



p = ^^|v+)<V+U®|o)<o|b + |^-)<V-U®|i)(i|b 
+|o)(ou ® liM<V>-lu + |i><iU ® I^+)^+IbJ . 

(3) 

where \ip±) = ^ (|0) ± |1)) and each qubit A and B have 
four nonorthogonal states. The two-qubit separable state 
([3|) show a product of quantum states that can not have 
any equivalent classical system and also present quantum 
correlation. Quantum discord quantifies all qauntum cor- 
relations including entanglement between A and B. 

Note that the mutual information T(A; B) is the 
sum of the quantum discord and the classical corre- 
lation CC(A;B) defined in ref. 14 , 1(A;B) = 6{A : 
B) + CC(A;B). In other words classical correlations 
CC(A;B) are equal to the mutual information J(A;B). 

Using quantum discord to discernate the quantum 
from the classical correlations we can study the behaviour 
of such different correlations in quantum phase transi- 
tions. In the next sections we focus our attention on 
the quantum phase transitions of the Ising and the XXZ 
spin- 1/2 chains and we study them using the quantum 
discord. 



III. QUANTUM DISCORD AND QUANTUM 
PHASE TRANSITION IN THE ISING CHAIN 

The Hamiltonian of the transverse one-dimensionnal 
Ising model reads 

N 

Hi stng = ~y] (Aof of +1 + erf ) , 

i=l 

with the boundary condition afj — af and a a with 
a = {x,y,z} are the Pauli matrices and cr° = I. Wc 
denotes by \g) = \ |) and |e) = | J.) the spin up and 
spin down states. For A = all spins are pointing in 
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the z direction while for A — > oo they point in the x di- 
rection. In the thermodynamic limit N — > oo the Ising 
spin chain undergoes a quantum phase transition at the 
critical point A c = 1. The correlation length diverges at 
this point^i. 

The quantum discord is worked out from the joint state 
of the two spins at the lattice sites i and j. The informa- 
tion on the joint state is contained in the two-site density 
matrix p^ which is derived from the following operator 
expansion 



1 3 

Pij = Tr [P] = T Qa 



0<?i 



a, /3=0 



where the coefficients Q a p of the expansion are related to 
the spin-spin correlation functions through the relation 



e aj£ 



Tr 



a 



o"; c^- 



Owing to the symmetry of the Hamitlonian most of the 
coefficients Q a /3 are equal to zero. Translation invariance 
require that the density matrix verifies pij — and 
is independent of the position i. The reflexion symme- 
try leads to p^ — pji, the Hamiltonian being real the 
density matrix verifies p*j — pij and the global rota- 
tion symmetry implies that the density matrix commutes 
with the operator of(T|. Combining all symmetry con- 
straints the density matrix expressed in the natural basis 
i\gg), \ge), \eg), \ee)} reduces tc£ 



Pij 




(4) 



with u± 



= i ± 



K 

4 2 

and y - 



KK) 

f — v 2 -) w z 

Kf£MfM) 
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i-KO 



The magnetiza- 



tion of the spin- 1/2 Ising chain is given byi^ 

(1 + A cos</>) 



and W0 = \J (A sin 4>) 2 + (1 + A cos </>) 2 is the energy spec- 
trum of the transverse Ising chain. The spin-spin correla- 
tions functions are related to the determinant of Toeplitz 



matrices 
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FIG. 1: Concurrence for nearest neighbour spins in the Ising 
spin 5=1/2 chain. The concurrence is maximal close to the 
critical coupling Ac = 1. 



and 
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'i+rl 



G r G- r where 



G K = — 

7T 



d<j> cos(<f>K) 



(1 + A cos < 



W0 



A C ax ■ ix ^ sin( ? !> 
/ a<psin(0Kj . 

71" JO W 

The quantum phase transition can be signaled by mea- 
sure of entanglement. Indeed the nearest and next- 
nearest neighbours entanglement reach their maximum 
at the critical point A^i A good measure of the entan- 
glement is provided by the concurrence C. The concur- 
rence of two spins may be computed from the joint state 
Pij through the formula C = max {0, 71 — 72 — 73 — 74} 
where the 7^ are the eigenvalues in decreasing order 
of the matrix R = J pijpfj 8 ' 9 - The matrix pfj is re- 
lated to the transpose of the two-site density matrix by 
pTj = {<J y ® <J y ) Pij [p v ®o v ). Hence the entanglement 
between two spins at site i and j is given by 

Cij = 2max{0, |x| - \y\ - w} . 

Figures [T] and [5] represent respectively the concurrence 
for the nearest-neighbour spins Cj^+i and next-nearest 
neighbour spins C^i+2. It shows that the maximum of 
entanglement is reached close to A c = 1 . Only the deriva- 
tives of the concurrence show singularities at the critical 
points' 5 signaling the presence of the quantum phase tran- 
sition. 

Entanglement increases close to the quantum phase 
transition. However entanglement is only one kind of 
quantum correlation. It is then legitimate to ask : What 
is the behaviour of the total amount of quantum cor- 
relations close to the critical point of the QPT ? The 
answer is provided by the quantum discord. The details 
of the derivation of the quantum discord are given in ap- 
pendix [X] Figure [3] represents the quantum discord for 
the Ising 5 = 1/2 chain for nearest-neighbour f^t+i and 
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FIG. 2: Concurrence for next-nearest neighbour spins in the 
Ising spin chain. The concurrence is maximal at the critical 
coupling A c = 1. 




FIG. 3: Quantum discord 8 for the nearest-neighbour (in 
dashed line) and next-nearest-neighbour (in full line) spin 
sites of the one-dimensionnal Ising model. The quantum dis- 
cord is minimized for <f) — and 9 varying progressively from 
zero for A < 1 to n/4 for A > 1, see appendix 1X1 for the details 
on the derivation. 



next-nearest-neighbour . It shows that the quan- 
tum correlations increase and are maximal in a region 
close to the critical point A c . 

More informations on the location and the order of the 
quantum phase transition can be obtained by consider- 
ation of the derivatives of the quantum discord with re- 
spects to the coupling parameter A. Indeed there is a re- 
lation between the nonanalyticity in energy and quantum 
discord that leads to the characterization of the quantum 
phase transition^. Consider the energy of the spin sys- 
tem in terms of the reduced density matrix of two spins 
at positions i and j, the energy reads 



^(p«) = E Tr [^. 



pij 



(5) 



FIG. 4: First derivatives of the quantum discord with re- 
spects to the coupling parameter A for the nearest-neighbour 
(in dashed line) and next-nearest-neighbour (in full line) spin 
sites of the one-dimensionnal Ising model. 



position i and j. The sum of Hij over the positions i,j is 
equal to the Hamiltonian of the spin system, H^ = 
H. For the one-dimensionnal Ising model the reduced 

Hamiltonian reads H*j in9 = — yXaf a™ + of / 'N 

where TV is the number of spins and 5ij is the Kronecker 
delta. The order as well as the location of a quantum 
phase transition can be characterized by the nonanalitic- 
ity of the energy. If the first derivative of the ground state 
energy presents a finite discontinuity then the quantum 
phase transition is of the first order. However if the first 
derivative is continuous while the second derivative shows 
discontinuity or divergence then the quantum phase tran- 
sition is of the second order. Moreover the derivatives of 
the energy are related to the derivatives of the reduced 
density matrix and one can show thati£ 




where Hij is the reduced Hamiltonian of the two spins at 



dE( Pij ) 
dX 

d 2 E{ Plj ) 
dX 2 



Note that the derivatives of the reduced Hamiltonian are 
continuous with respects to the coupling parameter A. 
Hence a discontinuity in the first derivative dE g f ^ 3 ^ of 
the energy at the critical point implies a discontinuity at 
least of one of the reduced density matrices pij . Similarly 
a discontinuity or a singularity in the second derivatives 

of the energy d requires the divergence of at least 

one of the derivatives at the critical points. It be- 
comes now evident that the quantum phase transition 
can be characterized by an analysis of the derivatives 
of the quantum discord with respects to the coupling pa- 
rameter A. The quantum discord being dependent on the 
reduced density matrix one deduces that i) a discontinu- 
ity in the quantum discord implies a discontinuity in /Oy 
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FIG. 5: Classical correlation CC for the nearest-neighbour 
(in dashed line) and next-nearest-neighbour (in full line) spin 
sites of the one-dimensionnal Ising model. 
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FIG. 6: Mutual information I(A; B) for the nearest-neighbour 
(in dashed line) and next-nearest-neighbour (in full line) spin 
sites of the one-dimensionnal Ising model. 



(hence the first derivative of the energy is discontinuous) 
and a first order QPT, ii) a singularity in the derivative 
of the quantum discord implies a discontinuity or a di- 
vergence of (the second derivative of the energy is 
discontinuous) and the QPT is of the second order. 

Figure [4] shows that for nearest neighbours ddi^ + i/dX 
presents a discontinuity at the critical point and for the 
next-nearest neighbours 881^+2/ d\ shows a singularity. 
The fact that the quantum discord is continuous while 
its derivatives are discontinuous or singular at A = A c 
indicates that the precise location of the critical point is 
at A c = 1 and that the QPT is of the second order. 

For A > 1 and approching the critical point from the 
right, the classical correlations decrease in opposite to 
the quantum discord which increases as depicted by the 
figures [3] and For A < 1 both classical and quantum 
correlations are decreasing when the coupling parame- 
ter A decreases. Classical correlations monotonically in- 
crease with respects to A while quantum discord show 
that the amplitude of the quantum correlations increase 
close to the critical point A c . From the mutual informa- 
tion we infer that the spins are uncorrelated for A — > as 
depicted in figure [5] and can be confirmed by the Hamil- 
tonian which reduces to Hi S i ng , a=o = — 52i=i a j- Close 
to the critical point the amount of quantum correlations 
is maximal as supported by the concurrence. 



IV. QUANTUM DISCORD AND QUANTUM 
PHASE TRANSITION IN THE XXZ CHAIN 



(ij). At the critical point A c = 1 the XXZ spin chain 
undergoes a quantum phase transition between an XY 
phase for — 1 < A < 1 and an Ising antiferromagnetic 
phase for A > 1. 

The density matrix of two nearest-neighbour spin sites 
is given by ([4]). However owing to the spin-flip symme- 
try the magnetization is equal to zero, (a z ) = and 
the XXZ Hamiltonian is also symmetric with respects 
to the rotation in the xy plan which leads to the equal- 
ity between the spin-spin correlation function in the xy 
plan (af<j'j) — (enVJ). Gathering all the symmetry con- 
straints of the XXZ Hamiltonian the joint density matrix 
Pij reduces to 



Pij 



( u 

w x 

x w 

V u 



with u± =u=\ + <CT ' 4 CTj) , w 
V = 0. 

The spin-spin correlation functions between nearest- 
neighbour spin sites for —1 < A < 1 are given b y 18 i 19 ' 20 



COS 7T$ 



dx 



; cosh : 



s'mhx cosh 2 ($2;) 

dx sinh ((1 — $)x) 



7r sinh(7r<£>) J_ co sinhx cosh $2; 
where A = cos (7r$), and 



The hamiltonian of the one-dimensionnal XXZ spin- 
1/2 chain reads 

Hxxz = ( S ? S i + S ? S j + AS i S l) > 
where the sum runs over the nearest-neighbour bonds 



2 . 2 f°° dx x cosh x 
(OiOi+i) = I--2 / Z^TZZZZz 



n"° J -co sinh x cosh (5>x) 
2cot(7r$) f°° dx sinh ((1 - $)x) 



sinh : 



cosh $a 



And for A > 1 the correlation function are given 

b 3 r 21 i 22 i 23 
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FIG. 7: Concurrence for nearest- neighbour spin sites in the 
XXZ spin S — 1/2 one-dimensionnal model. 
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FIG. 8: Quantum discord for nearest neighbour spins in the 
XXZ spin S = 1/2 one-dimensionnal model. The quantum 
discord is minimized for = and 6 varying progressively 
from zero for A < 1 to 7r/4 for A > 1, see appendix |A"1 for the 
details on the derivation. 



(« + l> = / — ^ 
J-oo+i/2 Sinfl(7TXj 

/ x , cot(^a;)\ 
x — — r cosh v — , 

V sin {(j>x) sinhy / 

with A = cosh v, and 

<^< + i> = 1+2 r +i/2 — 

J-oc+i/2 smh(Trx) 

x ( cot(^x) coth(V) J ) ■ 

V sin {vx) J 

The concurrence for the XXZ spin- 1/2 chain is given 
by C = 2max{0, \x\ — \u\}. Figure [7] represents the de- 
pendence of the concurrence with respects to the control 
parameter A and shows that entanglement between the 
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FIG. 9: Dependence of the classical correlation CC(A; B) on 
the parameter A for the nearest-neighbour spin sites in the 
XXZ spin 5 = 1/2 chain. 



nearest-neighbour spin sites is maximal at the critical 
point A c = l£. 

In parallel the quantum discord presents the maximum 
of the quantum correlations at the critical point. The 
derivation of the quantum discord for the XXZ model 
follow the same scheme as for the Ising model and the de- 
tails of the derivation are provided in appendixfAl Figure 
H] depicts the dependence of the quantum discord on the 
control parameter A. The quantum discord is minimized 
with respects to the angles <fi and 8 that control the axes 
of the projector ITj. 

As depicted in figure [5] the quantum correlations are 
stronger close to the QPT. On the opposite the clas- 
sical correlations are weaker close to the critical point 
A c . Both classical and quantum correlations measured 
by quantum discord and the mutual information J(A\ B) 
behave accordingly to the behaviour of entanglement 
measured by the concurrence. 

The cusps observed in figures [S] and [5] for both the 
classical correlation and the quantum discord arise from 
a level-crossing between the ground-state energy and an 
excited-state energy. This scheme explains the cusp ob- 
served in the concurrence measuring the entanglement 
of nearest neighbour spin in (d > 2)-dimensionnal XXZ 
model 4 -. A similar mechanism explains also the fact that 
the concurrence is maximal at the critical points for the 
one-dimensionnal XXZ model 24 . Level crossing induces 
a qualitative change in the ground states at the critical 
point A c = 1 and leads to the abrupt behaviour of the 
quantum discord and classical correlation at the quantum 
phase transition. As depicted in figure[7]the concurrence 
is less sensitive to the level-crossing than the quantum 
discord and the classical correlations since it does not 
present an abrupt cusp at A c . 

Remark that the quantum discord for the next-nearest 
neighbour spins is not presented for the Id XXZ 

model. The reason is that no more informations can be 
gained by than already furnished by the quantum 
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discord for nearest neighbours (Sf^+i. Indeed the quan- 
tum discord #i,j+2 presents a similar dependence on the 
coupling parameter A as shown by <5j,i+i and its ampli- 
tude verifies (5^+1 ~ 25^+2 ■ 
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V. CONCLUSION 

Quantum phase transition of the one-dimensional spin- 
1/2 Ising and antiferromagnetic XXZ models have been 
studied using the quantum discor d 10 ' 11 . Quantum dis- 
cord is a measure of the quantum correlation given by the 
difference of two classicaly equivalent expressions of the 
mutual information. Using quantum discord and mutual 
information we showed that it is possible to discernate 
the quantum and the classical nature of the correlations 
between two elements of a quantum many-body system. 

We demonstrated that the amount of quantum correla- 
tions is larger close to the critical points for the quantum 
phase transitions of the one-dimensionnal XXZ and Ising 
models. For the XXZ model the amount of classical cor- 
relations decrease at the critical point while for the Ising 
model they are simply monotonous with respects to the 
coupling parameter A. The behaviours of the quantum 
and classical correlations measured by means of the quan- 
tum discord are in agreement with the behaviour of the 
entanglement computed by the concurrence for both spin 
models. 

Quantum discord is a good way to discernate the na- 
ture of the correlations between the components of a 
quantum system and it is a good qualitative indicator 
of the existence of quantum phase transition. However 
depending on the model studied the maximum of the 
quantum discord is not necessarily located at the critical 
points as shown for the Ising model. The quantum dis- 
cord can only be used for qualitative detection of quan- 
tum phase transitions. On the contrary the derivatives 
of the quantum discord provide precise informations on 
the location and on the order of the quantum phase tran- 



sitions. Quantum discor d > and concurrence^'ii^ are 
two useful tools to detect quantum phase transitions. De- 
spite the fact that implementation of concurrence is more 
easy/rapid than for quantum discord, concurrence mea- 
sures only the nonlocal quantumness of the correlations 
between two systems A and B. Quantum discord mea- 
sures the total amount of quantum correlations between 
A and B including entanglement. 

We showed that quantum discord agrees with the pre- 
dictions provided by the concurrence for the behaviour of 
the quantum correlations close to quantum phase tran- 
sitions for one-dimensionnal quantum systems. It would 
be interesting to confirm this observation in futur investi- 
gation of quantum systems with higher dimensionnality. 
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APPENDIX A: DERIVATION OF THE 
QUANTUM DISCORD FOR SPIN CHAINS 

The quantum discord for two systems reads 5{A : B) = 
minr nB } {h {A) - H (A, B) + H (A\ {itf})} for which the 
three von Neumann entropy are computed over respec- 
tively the reduced density matrix pa, the joint state Pa,B 
and the conditional density matrix Pa\yi b ■ 

The reduced density matrix pa is given by 

PA = T^BPA,B 

= (u + + w)\g)(g\ A + (u-+w)\e)(e\ A , 

which leads to the von Neumann entropy of the reduced 
density matrix 



H(A) 



-Tr^p^log p A 



-id 



(O)log 
(O)log 



f 



(1 



(Al) 



The joint density matrix is equal to the density matrix 
given in equation ([4J pa,B = pij & n d the von Neumann 
entropy is given yb the sum of the Shannon entropy over 
the eigenvalues of the density matrix and reads 

H(A,B) = -Tta,bPA,B log p A ,B 

= — (w + x) log (w + x) — (w — x) log (w — x) 

-^S e logS £ , (A2) 



e=± 



where 



i{( U+ + n_)±[(u + -n_) 2 + V] 1/2 }. 



The conditional density matrix Pa\ti b is given by pro- 
jection on the arbitrary basis 

\9)b = oos(0)|ji) B + e^sin(0)|j 2 ) B , 
\e) B = e- < *sin(fl)|ii)B-coB(e)|j 2 )B. 

The angles <f> and 9 control the projectors direction and 
the quantum discord is minimized with repsects to these 
angles. 

The conditional entropy is then given by 



with Hf 

\ja)(ja\BPA,E 



-- 2.4 ® \ja)(j a \B and p ja = Tr A ,B^A 



8 



For the projection onto the state \j a ) where a = 1,2 
one gets 



\j a ){ja\B ® < |<?}(<?U^ Q ,+ + |e)(eUX Ja ,_ 



(A3) 



3J reads V jcn ± = -L-^(X jat+ + X joi} _) ± 

[(^, + -x ia) _) 2 + 4^ a y ia ] 1/2 |. 

Finaly the conditional von Neumann entropy reads 



where for a = 1, 



= (u + cos 2 (6») + wsin 2 (6»)) , 
X h - = (ii_sin 2 (6») +wcos 2 {6)) , 
Fj, = (e-^a; + 2/e^)sm(6l)cos(e), 

and with = | (1 + (o" z ) cos(20)). For the projection 
onto the state \j2) the elements Xj lt ± are given by 



x h,+ = (u+sm 2 (9) +wcos 2 (9)) , 
X,-,.- = (u- cos 2 (d) + wsm 2 (9)) , 



33,- 

Y, 



-- - (er^x + ye 1 *) sin(0) cos(0), 

= \ (1 - (cr z )cos(26»)). The eigen- 
values &j a ,± of the conditional density matrix 



and with 



Pj2 



H(A\{uf}) = ^ - Pj Tt A p A \ n n]ogp A \ n 3 

3 = 1,2 

e=± 3=1,2 



(A4) 



The quantum discord is minimized with respects to 
the angles cf> and 9 that control the direction of the 
projector ITj in the Hilbert space of the sub-system B. 
The minimization must be achieved for each ensemble 
of values of the magnetization and the spin-spin correla- 
tions function which defines the joint density matrix. It 
can be shown numericaly that for the one-dimensionnal 
Ising and antiferromagnetic XXZ models the set of an- 
gles {4>, 0} that minimize the quantum discord belong to 
the ranges 4> = 0, 9 € [0, 7r/4]. 
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